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Abstract 

Monte Carlo calculations of the irreducible particle-particle interaction on 
a two- leg Hubbard ladder doped near half- filling are reported. As the tem- 
perature is lowered, this interaction develops structure in momentum space 
similar to the magnetic susceptibility reflects the development of 

strong short-range antiferromagnetic correlations. Using this interaction, the 
eigenfunction of the leading singlet pair eigenvalue is found to have d^2_y2 
like symmetry. The single-particle spectral weight is also shown to peak near 
(vr, 0) and (0, vr) when the ratio of the inter- to intra-chain hopping t±/t ~ 1.5, 
leading to an increased tendency for pairing. 
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Numerical calculations, renormalization-group bosonization studies as well as strong- 
coupling treatments find that half- filled t-J or Hubbard two- leg ladders have a spin gapped 
ground state with short range antiferromagnetic correlations Furthermore, these vari- 

ous techniques all find that when holes are doped into the ladder, (i^.2_j^2-like pairing correla- 
tions develop. This being the case, one would like to understand the nature of the effective 
interaction that gives rise to the d^-i-y-i pairing correlations for this system. Here, in order 
to address this question, we present Monte Carlo results for the irreducible particle-particle 
interaction Fj. In addition, we use Fj and the single-particle Green's function obtained 
from the Monte Carlo calculations to solve the Bethe-Salpeter equation in the singlet pair- 
ing channel. We find that the eigenfunction with the leading eigenvalue has (ia.2_y2-like 
symmetry. 

The Hubbard model Hamiltonian for a two-leg ladder has the form 

i,X,s i,s iX 

Here t is the intra-chain one electron hopping, t± the inter-chain hopping and U the on-site 
Coulomb interaction. The operators cj^^ and create and destroy electrons of spin s on 
site i of the A*'^ leg, respectively, and n^s = c|as'^*As is the occupation number for spin s on 
site i of the A*'^ leg. 

Using Monte Carlo techniques, we have calculated the finite temperature two-particle 
Green's function 

G2(a;4,X3,X2,xi) = -{T c^{xi)ci{x3)cl{x2)c\{xi)). (2) 

Here cj(xj) with Xi = (x.j,rj) creates an electron of spin s at site Xj and imaginary time Xj 
and T is the usual r-ordering operator. Then, as previously discussed 0], one can take the 
Fourier transform of both the space and imaginary time variables and obtain G2{p', k', k,p) 
with p' = {p',iuJn'), etc. This two particle Green's function can be expressed in terms 
of the exact single-particle propagator Gs{p,iuJn) and the reducible particle-particle vertex 
T{p',k',k,p) 
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6*2 (p, k', k,p) = -5p^p> 5k,k' Gi{k)G^{p) (3) 
+^^'+k',p+k G^{p')G^{k')T{p', k\ k,p)G^{k)G^{p). 

Then from the Monte Carlo data for G and G2, one can determine r{p' , k' , k,p). Finally, 
because the effective pairing interaction corresponds to the irreducible particle-particle in- 
teraction Fi in the zero energy and momentum center of mass channel, we have inverted the 
fully dressed t-matrix equation to find Fj in terms of F and G. Setting k = —p and k' = —p', 
the fully dressed t-matrix equation becomes 

r{p'\p) = r,{p'\p) -LY^r{p'\k)G^^{-k)G^{k)r,{k\p), (4) 

k 

which is solved to find Fi(p'|p). The effective pairing interaction in the singlet channel is 

Vip'-p) = l{Up'\p) + U-p'\p)). (5) 

We have carried out this calculation at different temperatures for various values of the 
hopping anisotropy t±/t, interaction strength U/t and filling (n) = {rii^ + riii). Here, we 
will show results for tj_/t = 1.5, U/t = 4, {n) = 0.875 and a 2 x 16 lattice. We have 
chosen this set of parameters, because the numerical density matrix renormalization group 
(DMRG) calculations find that for U/t = 4 and (n) = 0.875, the d^-^-y-i pairing correlations 
are strongest when t±/t ~ 1.5 p. In addition, for these parameters we have good control 
of the maximum entropy analytic continuation of the Monte Carlo data, which is necessary 
to obtain the single-particle spectral weight y4(p, cu). 

In Fig. |l| we plot the effective pairing interaction V versus = v'x ~ Vx for = 
p'y — Py = vr. Here we have set uJn = '^n< = ttT corresponding to Um = energy transfer. 
The three curves correspond to temperatures T = l.Ot, 0.5t and 0.25t. One can see that as 
the temperature decreases, the effective pairing interaction becomes increasingly positive at 
large momentum transfer q — > (tt, tt). As shown in Figure |^, the magnetic susceptibility 

X(q) = ^ I' dr Y: e'^-^ (m-+,(r)m+(0)) (6) 
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also develops structure at large momentum transfers over this same temperature region. 
Thus it is clear that the effective pairing interaction is associated with the development of 
short-range antiferromagnetic correlations. 

It is also of interest to study the Bethe-Salpeter equation for the singlet particle-particle 
channel 

T 

A„0a(p,iu;n) = -— ri(p,it^„|p',iu;„/)|G(p',iu;„/)p0o(p',iu;„/). (7) 

Fig. ^ shows the eigenfunction 0(p, iuOn) of the leading eigenvalue versus p for uOn = ttT. We 
observe that 0(p, iiiT) has dx2-y2-\ike momentum structure in the sense that it has opposite 
signs and is largest near (tt, 0) on the bonding band and near (0, vr) on the antibonding band. 
This is associated with the structure of the irreducible interaction Fi and the single-particle 
spectral weight A{p,Lj) which we will study below 0. 

The temperature dependence of the leading eigenvalue Ai is plotted in Fig. H(a). As 
previously discussed, when the temperature is lowered, short-range antiferromagnetic corre- 
lations develop and the effective pairing interaction increases at large momentum transfer. 
This leads to an increase in Ai as shown in Fig. §(a). We believe that Ai will approach unity 
at low temperatures where power-law ci^^z-ya-like pairing correlations have been shown to 
exist using density matrix renormalization group techniques 0. In Fig. §(b), we show the 
dependence of the leading eigenvalue Ai on the hopping anisotropy tj_/t. According to the 
DMRG calculations 0, for U/t = 4 and (n) = 0.875 the pairing correlations in the ground 
state are strongest when tj_/t ~ 1.5. In Fig. §(b), we do not observe a strong dependence 
of Ai on t±/t because of the thermal smearing effects at T = 0.25t. 

In addition to the irreducible interaction vertex Fi, the single-particle Green's function 
G(p, iun) is also important in determining the structure of the leading eigenfunction of 
the Bethe-Salpeter equation. Using a numerical maximum entropy procedure [0], we have 
calculated the single-particle spectral weight 

A(p,uj) = --lmG{p,u}). (8) 

TT 
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This is plotted in Fig. ^ for t±/t = 1.5 as a function of u for different p. Here, the sohd 
curves are for the bonding band {py = 0) and the dotted curves are for the antibonding band 
{py = it). We see that for this value of t±/t, the bonding band has spectral weight near 
the Fermi level for p ~ (tt, 0), and the antibonding band has spectral weight near the Fermi 
level for p ~ (0,7r). Hence, these Fermi points can be connected by scatterings involving 
q = (tt, tt) momentum transfer. Since Fjs is large and repulsive for q ~ (vr,7r), the leading 
eigenfunction of the Bethe-Salpeter equation has opposite signs for p near p = (vr, 0) and 
(0, vr), as seen in Fig. ||. 

These calculations provide further insight into the structure of the effective pairing in- 
teraction and the single-particle spectral weight which lead to the pairing correlations in the 
two-leg Hubbard ladder. Specifically, the momentum structure of the effective interaction 
V{q) clearly reflects the existence of short-range antiferromagnetic correlations as the cause 
of the increasing positive strength of V^(q) at large momentum transfer. Secondly, the en- 
hanced spectral weight in the bonding band near (vr, 0) and the antibonding band near (0, vr) 
are reminiscent of a similar effect observed in the two-dimensional Hubbard model near half- 
filling 1^,^] and in the ARPES of the cuprates |]TD[ . This enhanced low lying spectral weight 
associated with the renormalized quasiparticles is such that pair scattering processes with 
momentum transfers near (vr, vr) have increased phase space. These two features appear to 
play an important role in the development of pairing on the two-leg ladder, just as they do 
for the two-dimensional Hubbard model. 
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FIGURES 




FIG. 1. Momentum dependence of the effective interaction V^(q) for U = 4i, (n) = 0.875 and 
t± = IM. Here I^(q) is measured in units of t, qy = tt and V{q) is plotted as a function of Qx- 
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FIG. 2. Momentum dependence of the magnetic susceptibihty x(q) for U = At, (n) = 0.875 
and t± = IM. Here Qy = and x(q) is plotted as a function of Qx- 
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FIG. 3. Momentum dependence of the d^2_y2 eigenfunction 0(p, iTrT). These results are for 
T = 0.25t, U = At, {n) = 0.875 and t^. = IM. 
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FIG. 5. Single-particle spectral weight A{p,u}) versus lo for t±/t = 1.5, T = 0.25t, U/t = A 
and (n) = 0.875. The solid curves denote the results for the bonding band {py = 0) and the dotted 
curves denote the results for the antibonding band {py = tt). 
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